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SUMMARY 


This paper quantitatively examines the precision of the model of pattern motion 
analysis put forth by Adelson and Movshon (1982) who proposed that humans 
determine the direction of a moving plaid (the sum of two sinusoidal gratings of 
different orientations) in two steps. The velocities of the grating components are 
first estimated, then combined using the intersection of constraints to determine the 
velocity of the plaid as a whole. Under the additional assumption that the noise 
sources for the component velocities are independent, an approximate expression 
can be derived for the precision in plaid direction as a function of the precision in 
the speed and direction of the components. Monte Carlo simulations verify that the 
expression is valid to within 5% over the natural range of the parameters. The 
expression is then used to predict human performance based on available estimates 
of human precision in the judgment of single-component speed. Human 
performance is predicted to deteriorate by a factor of 3 as half the angle between the 
wavefronts (0) decreases from 60° to 30°, but actual performance does not. The 
mean direction discrimination for three human observers was 4.3° + 0.9° (SD) for 
0 = 60° and 5.9° ± 1.2° for 0 = 30°. This discrepancy can be resolved in two ways. If 
the noises in the internal representations of the component speeds are smaller than 
the available estimates or if these noises are not independent, then the 
psychophysical results are consistent with the Adelson-Movshon hypothesis. 

INTRODUCTION 

Adelson and Movshon (1982) proposed that the human visual system 
determines the direction of a moving pattern using a two-step procedure whereby 
the velocities of oriented components within the pattern are first estimated then at a 
later stage are combined to calculate the motion of the pattern as a whole. This 
hypothesis was originally formulated to explain their psychophysical findings that, 
in order for a plaid (the sum of two sinusoidal gratings of different orientations) to 
be seen to move coherently as a whole, the component gratings must be similar in 
both spatial frequency and contrast. They concluded from these psychophysical 
results that the human visual system analyzes plaid motion by first decomposing it 
into the motion of the grating components (Fig. 1). They suggested that this 
decomposition is the natural consequence of having orientation and spatial- 
frequency tuned motion sensors at the front end of the system (for a review, 
DeValois and DeValois, 1980). They also proposed that, at a second stage, the 
component velocities are recombined using the intersection of perpendicular 
constraints (Fennema and Thompson, 1979) to yield a measure of the motion of the 
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plaid as a whole. The lack of coherence for gratings of widely differing spatial 
frequencies was explained by assuming that the second stage of the process is only 
possible among sensors with similar spatial-frequency tuning. 

They found support for their hypothesis in the discovery of two types of motion- 
sensitive neurons in the monkey visual cortex: one sensitive to component motion 
and one, at a higher level within the cortex, sensitive to the motion of the plaid as a 
whole (Movshon et al., 1986). Furthermore, a recent study has found that speed 
discrimination for moving plaids is consistent with the component-driven 
approach that Adelson and Movshon proposed (Welch, 1988). 

Recently, however, there have been a number of challenges to the Adelson- 
Movshon hypothesis. Two studies have shown that when the gratings differ in 
contrast that the perceived direction of motion deviates significantly from that 
predicted by the intersection of constraints (Kooi, DeValois, and Wyman, 1988; 
Stone, Mulligan, and Watson, 1988). This problem, however, can be overcome if 
the perceived speed of the grating components is a function of contrast (Thompson, 
1982) and perceived grating speed rather that actual grating speed is the input to the 
second stage of processing. In addition, even for plaids with gratings of equal 
contrast, certain combinations of grating speeds and orientations yield perceived 
directions of plaid motion significantly different from that predicted by the 
intersections of constraints (Ferrera and Wilson, 1988). 

This study examines quantitatively the theoretical consequences of the Adelson- 
Movshon hypothesis on the discrimination of the direction of a moving plaid. An 
expression for the uncertainty in plaid direction as a function of the uncertainty in 
component velocity is derived assuming that the channels carrying the component 
information have exclusively independent sources of noise. Human performance 
shows that the direction of plaid motion is often determined with greater precision 
than is expected from the precision in independent single-component judgments 
within the context of the Adelson-Movshon model. These results indicate that, if 
the Adelson-Movshon hypothesis is correct, either the available measures of speed 
discrimination (McKee, Silverman, and Nakayama, 1986) are poor estimates of the 
internal noise associated with single-component speed or the channels carrying the 
component speed information share at least some of their noise sources. 

The author would like to thank Al Ahumada, Jeff Mulligan, John Perrone, and 
Andrew Watson for insightful conversations and comments on an earlier draft. 
This work was partially supported by an NRC postdoctoral associateship. 
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Figure 1. - Component-oriented reconstruction of pattern velocity. Adelson and Movshon (1982) proposed 
that plaid velocity was reconstructed from the velocity of the two grating components. The measurement 
of component speed is known to have some inherent noise associated with it. This noise, shown as the 
perturbations si and s2, would be expected to propagate through the system and generate noise in the 
determination of plaid velocity. Noise in the determination of grating orientation has been neglected. 
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THEORY 


Direction of Plaid Motion as a Function of Component Motion 

The direction of motion of a plaid is a function of the relative orientation of the 
two grating components and their respective speeds. In the following analysis, I 
will determine that function. By definition, each grating component is oriented 
symmetrically with respect to the bisecting axis, i.e., each grating is offset equally 
with respect to the axis that symmetrically divides the plaid in half. The following 
calculation of direction is made with respect to this bisecting axis which is assumed, 
for simplicity, to be vertical. Figure 2 illustrates a moving plaid in velocity space 
with the horizontal axis indicating horizontal speed and the vertical axis indicating 

vertical speed. The two vectors (v^ and ) represent the velocities of the two 

grating components while the vector \T represents the velocity of the plaid as a 

whole. The two components move symmetrically, 8 degrees from vertical, at 
speeds of Vj and V 2 , respectively, while the plaid as a whole moves at an angle a 
with respect to vertical, at a speed of Vp. 

From figure 2, the intersection of perpendicular constraints (Fennema and 
Thompson, 1979) can be used to relate the motion of the plaid as a whole to the 
motion of the individual grating components. Specifically, plaid motion can be 
decomposed into component motion simply by calculating the orthogonal 

projections of \T onto the 6 and -6 directions, yielding 


and 


Vl=Vp cos (6 - a) 

(la) 

V 2 = Vp cos ( 0 + a) 

(lb) 
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Figure 2. - Diagram illustrating plaid motion in velocity space. V] and V 2 are the speeds of the 
component gratings perpendicular to their respective wavefronts (in the directions 6 and -6, 
respectively, with respect to the bisecting axis). Vp is the speed of the plaid as a whole; a is the 
direction of plaid motion with respect to the bisecting axis. The bisecting axis in this case is the 
vertical axis Vy. 

Combining these two equations to eliminate Vp, one can show that the direction of 
plaid motion with respect to the bisecting axis (a) is 


a = arctan 


cotan 9 


v,-v 2 

[V.+vJ 


( 2 ) 


A complete derivation can be found in the appendix. Note that if Vj = V 2 , plaid 
motion is along the bisecting axis and a is zero. 


Noise in the Direction of Plaid Motion Induced by Noise in Component Motion 

To determine the effect of noise in component velocity on the direction of plaid 
motion, I first determine the effect of small perturbations in component speed (sj 

and S 2 ) and component direction (5j and 5^) on the direction of plaid motion (a). 
Just as above, plaid motion can be decomposed into component motion as follows: 
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and 


Vi + Sj = Vp cos ( 0 + 5i - a) 
V 2 + s 2 = Vp cos (0 + S 2 + a) 


(3a) 


(3b) 


where equations (3) are derived from equations ( 1 ) with s\, sj, 8\, and ^ added as 
small perturbations in Vj, V 2 , and 6 . Remember that each grating is oriented 

symmetrically with respect to the bisecting axis so 6 is the same for both gratings. 
However, the noise in the direction of motion of the two gratings is assumed to be 
independent (5| is independent of 82 ) as is the noise in the speeds of the two 
components (S| is independent of S 2 ). 

Combining the two equations to eliminate V p , one can show that the direction of 
plaid motion (a) is 


a 


F(sj , S 2 , <5i , 82 ) = arctan 


(Vj + Si ) cos (0 + 82 . ) - (V 2 + S 2 ) cos (0 + <5i ) 
(V! + Si ) sin (0 + 81 ) + (V 2 + S 2 ) sin (0 + Si ) 


(4) 


Note that equation (4) reverts to equation (2), if si = S 2 = 5i = 82 = 0. 


The problem of determining the effect of component noise on plaid direction 
becomes that of determining the standard deviation of a (<r) as a function of the 
standard deviations of si, S 2 , and ^ ( <r Sl , ct S2 / o - 5 1 , o§ 2 ). This task can be 
made simpler if the nonlinear function F in equation (4) is approximated by the 
linear portion of its Taylor series: 


9F 3F 3F 9F 

a = F(si,s 2 , 8 i, 8 z) =F(0,0 / 0,0) + — Si + — S 2 + — — $ (5) 

3si ds 2 d5z 

with F(0,0,0,0) being the mean directional bias as defined in equation (2) and the 
remaining terms being perturbations on that mean. Since F is now a linear 
function, if Sj, s 2 , 8^, and ^ are independent random variables with Gaussian, 

zero-mean, distributions, then the standard deviation of a is given by the 
following formula: 


6 



a 


( 6 ) 



aF 


\2 


dsi 


Os 2 + 


fMSP 


Therefore, to determine o , one must first take the partial derivatives of F (eq. 4) 
with respect to sj, S2, and , evaluated at (0, 0, 0, 0): 


dF 2 V 2 tan 9 

dsi (V 1 +V 2 ) 2 tan 2 G + (V^V 2 f 

dF _ -2 Vi tan B 

as 2 (V 1+ V 2 f tan 2 0 + (Vt-V,) 2 
0F _ V 2 (V 2 -Vi) + V 2 (Vi + V 2 )tan 2 0 

(Vi + V 2 ) 2 tan 2 0 + (V r V 2 ) 2 
aF = Vi(V 2 -Vi) - Vi (Vi + V 2 ) tan 2 0 
dS 2 (Vi +V 2 ) 2 tan 2 Q + (Vi V 2 ) 2 


(7) 

( 8 ) 

(9) 

(10) 


Combining equations (7), (8), (9), (10) with equation (6) yields the general formula 
for the noise in plaid motion as a function of the noise in component motion. A 
complete derivation is given in the appendix. 

Let us now look at the special case in which the plaid is moving straight up (i.e., 
Vl=V 2 = V) and the noise amplitudes for each component are equal (i.e., 

<7 S = cr Sl = cr^ and 05 = crg 1 = <r§ 2 ), the general formula reduces to 


a _ V g s 2 + V 2 tan 2 8 cl (11) 

V tan0 V2 

Furthermore, if one assumes that the noise in V is proportional to V (McKee, 1981) 
and that A is the normalized standard deviation, i.e.. 


<t s = AV, 


( 12 ) 
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then equation 11 further reduces to 


C 



+ tan 2 9 Qg 


tan 6 V2 


(13) 


One now has an equation relating the standard deviation in plaid direction to the 
standard deviations in component speed and direction. 


Testing the Linear Approximation 

Equation (13) was derived using a linear approximation (i.e., eq. 5). To verify the 
validity of this approximation, I compare results generated by equation (13) with 
those generated by Monte Carlo simulations of equation (4). The results are shown 
in table 1. The third column is a given by Monte Carlo simulations of equation (4) 
using 10,000 iterations setting the perturbations si, S 2 , and ^ to random 
variables with independent, zero-mean, Gaussian distributions. The fourth column 
is a given by equation (13). Note that, although all values in this table are given in 
degrees, angular values for a§ in equations (11) and (13) must be in radians and 

equations (11), (13), and (14) yield values of a in radians. 


The agreement between third and fourth columns in table 1 clearly illustrates 
that the linear approximation was justified and that, for an extensive range of 6 , 

A, and o s , equation (13) provides a good estimate of the actual standard deviation 
of a. The first three rows confirm the fact implicit in equation (13) that, if the noise 
in V is proportional to V, then a is independent of V. The remaining rows show 
that changes in 0 , A, and a s cause significant changes in a, but that agreement 
between the third and fourth columns is within 5% and usually better. 
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TABLE 1. - VALIDATION OF THE LINEAR APPROXIMATION 


Conditions 

Speed (V), 

deg/9ec 

Monte Carlo, 

deg 

Equation 13, 

deg 

% error 

& 

II 

<2> 

o 

t—* 

II 

& 

O 

t-H 

II 

0.1 

2.46 

2.44 

0.8 


1 

2.45 

2.44 

0.4 


10 

2.48 

2.44 

1.6 

1 

II 

CD 

o v 
CM 

II 

s' 

o 

r— ' 

II 

1 

2.76 

2.73 

1.1 

X = 10%, <7g =10°, 0 = 60° 

1 

7.54 

7.45 

1.2 

X = 20%, <T§ = 1°, 0 = 60° 

1 

4.97 

4.73 

4.8 

X = 10%, <r§ = 1°, 0 = 45° 

1 

4.10 

4.11 

0.2 

X = 10%, = 1°, 0 = 30° 

1 

7.12 

7.05 

1.0 


Prediction of Human Direction Discrimination 

If we assume that the measured discrimination of single grating speed (McKee, 
Silverman, and Nakayama, 1986) and orientation (Caelli et al., 1983; Regan and 
Beverley, 1985) are limited by the same internal noise (A and <J S ) as the model in 
figure 1, then one can predict the precision in the direction of plaid motion from the 
precision in grating speed and direction using equation (13). To do so, however, 
requires knowing the precision of single-component judgments. 

The best human performance in speed discrimination is around 5% (McKee, 
1981). However, for a 1.5 cycle/deg sinusoidal grating moving at l°/sec, a more 
appropriate value is about 11% (Welch, personal communication; McKee, 

Silverman, and Nakayama, 1986). Since the above values of speed discrimination 
are measures of the semi-interquartile difference (half the range between 25 and 75% 
performance), A , the normalized standard deviation (see eq. 12), is 1.48 times larger 
or about 16% for the stimulus in the experiment described below. 

Direction discrimination for moving random dot patterns is about 1° or 2° 
(DeBruyn and Orban, 1988). For moving gratings, however, direction and 
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orientation discrimination are confounded. Because of this, it is difficult to 
determine how well the direction of moving gratings might be discriminated in the 
absence of orientation cues. However, orientation discrimination for a 1.5 cycle/ deg 
grating oriented 30-60° from vertical is probably 2° or higher (Caelli et al., 1983; 
Regan and Beverley, 1985). In any case, examination of table 1 (see second and 
fourth rows) shows that for values of A of 10% (or higher) that X dominates over 
values of a s near 1°, i.e., doubling Og produces only a 10% change in plaid 
direction discrimination. Therefore, the precise value of Og is not critical. For 
simplicity, the small contribution to the uncertainty in the direction of plaid motion 
from a g is neglected altogether and is omitted from the flow chart in figure 1. 

This yields the following simplified formula: 



Y2 


Equation (14) is an underestimate of the uncertainty in the direction of plaid 
motion because it assumes no noise in the internal representation of component 
direction and because it assumes that no noise is associated with the second stage of 
processing. Furthermore, X = 16% is a conservative estimate of fractional speed 
discrimination in that it was measured under ideal circumstances (high contrast, 
long duration, large number of cycles, and highly practiced subjectsXMcKee, 
Silverman, and Nakayama, 1986). Equation (14) also assumes that the noise sources 
associated with the internal representations of component speed for the two gratings 
are totally independent. If the noises in component speeds are correlated and 
proportional to speed (McKee, 1981), they will produce noise in plaid speed, but no 
noise in plaid direction. This can be seen graphically by noting that, if Vj and V 2 in 
figure 2 are incremented or decremented by the same percentage, there is no change 
in plaid direction. This can also be shown analytically by noting that, if V| and V 2 

are changed by the same proportion in equation (2), a remains unchanged. 
Therefore, if there is some correlation in the component noise, the predicted 
direction discrimination will be given by the following formula: 

a = M cotane (15) 

V2 

with U0, the fraction of the noise in the component speeds that is uncorrelated. 
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EXPERIMENT 

Methods 

To test the prediction of the previous section, human direction discrimination 
was measured with moving plaids. Three subjects viewed a plaid moving at 2°/s 
for 300 msec in some direction between 30° to the left and 30° to the right of straight 
up and were asked whether the plaid moved to the right or to the left of subjective 
vertical. Subjects were instructed to fixate a small cross which appeared in the 
center of the screen 500 msec before the stimulus. The fixation point was 
extinguished during the stimulus presentation. The plaid stimulus was composed 
of two 1.5 cycle/deg sinusoidal gratings oriented symmetrically with respect to the 
vertical axis with 6 equal to 60, 45, or 30°. The direction of plaid motion was 
changed by changing the relative speeds of the two grating components such that 
plaid speed (V p ) and orientation were always the same. The plaid image was 
displayed within a Gaussian window (spread = 0.25 image width) on a 19-in. 
monochrome monitor using an Adage RDS-300 graphics system. The stimulus was 
composed of superimposed halftoned gratings and was made to move using a 
lookup-table animation technique described in detail elsewhere (Mulligan and 
Stone, 1988). Subjects viewed the 5.3° by 5.3° image from a distance of 273 cm 
(thereby making the halftoning noise invisible) and used a chin rest to stabilize their 
heads. A simple updown staircase method (Cornsweet, 1962) was used to select the 
direction for a given trial. Two independent, randomly interleaved staircases were 
run for each value of 0. 

When human observers were asked whether the plaid was moving to the right 
or to the left of the vertical axis, their performance generated a sigmoidal curve in 
plots of "percent of trials where the plaid is seen to move to the left" versus "actual 
direction of plaid motion." When the plaid moved strongly to the left, observers 
always responded "left." When the plaid moved strongly to the right, observers 
never responded "left." In between, there was a range of directions for which 
observers had a varying degree of uncertainty: they responded "left" a certain 
percentage of the time. 

Human psychophysical data generated using this procedure was well fit by a 
cumulative Gaussian curve (the integral of a Gaussian). A simple way of 
quantifying human performance is therefore to determine the best fitting Gaussian 
(Finney, 1971). Then, one can use the mean as a measure of perceived vertical (the 
direction for which subjects responded "left" 50% of the time) and the standard 
deviation as a measure of direction discrimination (the range of directions over 
which there is uncertainty). The latter is plotted as o in figure 3. 
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Results 


Direction discrimination for moving plaids is often better than that predicted by 
equation (14), but is well predicted by equation (15) using the simple assumption 
that u 0 = sin© . Figure 3 shows a plot of o for the three subjects, from equation 
(14), and from equation (15) for three different values of 0. At 0 = 60°, human 
performance (mean a = 4.3°) is well predicted by both equation (14) (3.7°) and 
equation (15) (3.2°). However, all three subjects showed only a small increase in or 
(decrease in performance) with increasing 6 . This is as predicted by equation (15) 
but equation (14) erroneously predicts a threefold decrease in performance. So at 
0 = 30°, the mean human performance (5.9°) is well predicted by equation (15) 

(5.6°), but is on average nearly twice as good as predicted by equation (14) (11.2 °). 
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Figure 3. - Plot of the plaid direction discrimination of three subjects (1 and 2 were ignorant of the goals 
of this study) and two different predictions of a versus the angle between the component wavefronts 
(6). Psychometric curves were fit with a cumulative Gaussian whose standard deviation (a) is 
plotted here as a measure of direction discrimination. Each data point is the mean of three runs with 
each run consisting of 64 observations per psychometric curve. The standard deviation for each point 
was always less than 1.6° and averaged 0.85° and those for subject 3 are shown explicitly. The points 
labelled equation 14 and 15 were evaluated using X = 16%. 
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DISCUSSION 


Equation (15) can explain human direction discrimination of plaid motion if one 
assumes that the noises in grating speeds are correlated in relation to the angle 
between the gratings. Although there is no direct evidence for this type of 
correlation, it seems reasonable to assume that neurons whose preferred directions 
are close might share a significant amount of their input and therefore a significant 
amount of their noise. Assigning sin© as the fraction of noise that is uncorrelated 
was arbitrary. Further studies over a wider range of angles will be necessary to get a 
good estimate of uq. 

The predicted curve in figure 3 generated by equation (14) is both too high and 
the wrong shape. It underestimates human performance and is very sensitive to 
the relative orientations of the gratings within the plaid ( 6 ) whereas actual 
performance is not. However, an important caveat must therefore be kept in mind 
when interpreting this result: the predicted curve relies on the assumption that X 
is 16% or higher. This assumes that the noise in the internal estimate of speed used 
by the second stage is expressed uncorrupted as uncertainty in judgments of single 
grating speed (McKee, Silverman, and Nakayama, 1986). However, it is possible that 
the single-component judgments made by McKee and colleagues are noisier than 
the internal estimates of component speed that the second stage of the model in 
figure 1 has access to. In other words, X may be lower than 16%. Lower values of 
X could move the entire predicted curve down below the data. Flexibility in the 
choice of X can therefore make the predictions generated by equation (14) at least 
compatible with the psychophysical results although other noise sources at later 
stages would be necessary to flatten the curve. 
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CONCLUSION 


Human discrimination of the direction of a moving plaid can be quantitatively 
predicted by the model proposed by Adelson and Movshon (1982), although it 
constrains how the model might actually be implemented with the visual cortex. 
The direction discrimination data presented here can be explained in at least two 
ways. One possibility is that the uncertainty in component speed is lower than 
suspected (McKee, Silverman, and Nakayama, 1986) and not rate limiting. If so, 
direction discrimination would be limited by some noise source at a later stage of 
motion processing or in the decision process. Another possibility is that the 
uncertainty in component speed is the rate limiting noise, but performance is 
enhanced and made robust to changes in 6 by correlation in the noise between 
channels carrying information about the individual components. Both of these 
explanations assume that the model in figure 1 is essentially correct. Unfortunately, 
human direction discrimination for moving plaids does not resolve the issue of 
whether or not the Adelson-Movshon hypothesis is correct. 
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APPENDIX 


1. Proof of equation (2): 


From equations (1), eliminating Vp , 


cos(0 - a ) cos(0 + a ) 


then expanding the cosine terms, 


cos0 cosa + sin0 sina cos0 cos a -sin0 sina 


regroupmg. 


oos0 cosa 


(-1 l-\ = sin0 sina (--1 L 

\Vi V 2 / l Vi V 2 , 


tana =SQSfi. ( Yl _ Y? . | = cotan0 f Vl ~ V ^ 


sin0 - -1 — L 
V Vi V 2 ( 

which immediately leads to equation (2). 


Vi+V 2 


2. Proof of equation (4): 


Similarly, from equations (3), 


Vi + si 


V 2 + s 2 


cos(0 + 5i - a ) cos (0 + 5i + a ) 


Vi + si V 2 + s 2 

cos (0 + 5i) cosa +sin(0 + 5i) sina cos(0 + S2) cosa -sin(0 + 5 2 ) sina 


sin(0 + Si) sin(0 + 82) 1 cos(0 + 82) cos (0 + £ 1 ) 

sina + 1 = cosa 

* Vi + si V 2 + s 2 / V 2 + S2 Vi + si 


tana = 


(Vi + si) cos(0 + 82) - (V 2 + S 2 ) cos(0 + £ 1 ) 
(Vi + si) sin(0 + S 2 ) + (V 2 + s^ sin(0 + 5i) 


which immediately leads to equation (4). 


3. Derivation of equation (7) and (8): 


From the chain rule. * f 2 (x)f 

dx dx 


so from equation (4) 
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/ \2T"1 

1+ I (Vl + s l) cosjO + 62) - (V 2 + S2> cosjd + <5i) I 

\(Vi + si) sin(0 + 82) + (V 2 + S 2 > sin(0 + 8j)f 


and B = 

setting Sl = S2a g 1 = g 2 = 0 / jI= [1 + f ^ 1 ~ v 2 )cosg j 2 ] 1 r_2o^g_anfl V 2 ' 

dsi L m + V 2 ) sinfl / J [((V! + v 2 ) sine f_ 

3F 
3si 

and dividing through by cos 2 0 yields equation ( 7). 

By symmetry (i.e., exchanging the subscripts and multiplying by -1), one can generate equation (8). 

4. Derivation of equations (9) and (10): 

As above, = A B 

38 1 

with A unchanged and 

(V2 + S 2 )sin(fl + <5i)((Vi + si) sin(0 + fo) + (V2 + S2) sin(e + 0i))-(V2 + S2)cos(0 + 0i) ((Vi + si) cos(e + 6 2) - (V2 + S2) cos(e + 5i)) 

((Vi + si) sin(0 + 1S2) + (V2 + S2) sin(0 + Si)f 

setting s 1 = S2 = 5i = 52 = 0, = 

38 1 

9F _ [ sin 2 e V 2 (Vi+V 2 )-cos 2 e V 2 (Vi -V 2 ) 
a<Si [ (Vi+v 2 ) 2 sin 2 e + (Vi - V 2) 2 cos 2 e 

Dividing by cos 2 0 yields equation (9); exchanging subscripts and multiplying by -1 yields equation (10). 


-V2)cose \ 2 [ sin 2 e V2 (Vi + V2) - cos 2 e V2 (Vi - V2) 


(Vl+ V2)sine 


((Vi + V 2 ) sine f 


2 cose sine V2 

(V1+V2) 2 sin 2 e + (Vi -V 2 ) 2 cos 2 e 


cos (0 + 02)((Vi + »l) sin(0 + +JV 2 + S2) sin(0 + 0i))-sin(O + 02) ((Vi + si)cos(e + 02)-(V2 + S 2 ) eos(0 + gp) 

((Vi + si) sin(0 + 82) + (V2 + S2) sin(0 + Sj)f 
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